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Abstract 

Many authors claimed that a large initial inhomogeneity prevents the onset 
of inflation and therefore inflation takes place only if the scalar field is homo- 
geneous or appropriately chosen over the horizon scale. We show that those 
arguments do not apply to topological inflation. The core of a defect starts 
inflation even if it has much smaller size than the horizon and much larger 
gradient energy than the potential, as long as the vacuum expectation value 
is large enough 0.3mpi) and the core is not contracting initially. This is 
due to stability of false vacuum. 
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I. INTRODUCTION 



Inflation gives a natural solution of the horizon problem of the big-bang universe. How- 
ever, if inflation requires homogeneous initial conditions over the super-horizon scale, one 
cannot say that it is a solution of the horizon problem, though it reduces the problem by 
many orders of magnitude. To answer this naive question, many people have investigated 
how initial inhomogeneity affects the onset of inflation [1-7]. The most remarkable work was 
done by Goldwirth and Piran [1], who solved the full Einstein equations for a spherically 
symmetric spacetime. Their results are summarized as follows. New inflation is so sensitive 
to initial inhomogeneity that it requires homogeneity over a region of several horizon sizes 
is needed. Chaotic inflation is not so affected by initial inhomogeneity but requires a suf- 
ficiently high average value of the scalar field over a region of several horizon sizes. Later 
many people explored the same problem in different approaches [3-7] , which as a whole con- 
firmed the results obtained by Goldwirth and Piran. This suggests failure of the "ordinary" 
scenario of inflationary cosmology and therefore some alternatives were discussed [5-7], 

In this paper we investigate whether topological inflation [8] also requires initial homo- 
geneity over the horizon size. The simple model of topological inflation is described by the 
potential. 



However, the same type of inflation takes place in many plausible theories such as a nonmini- 
mally coupled massive scalar field [9], some models of supergravity [10], and some superstring 
inspired models [11,12]. The inflationary condition for the model (1.1) was investigated by 
numerical analysis [13]: r) ^ 0.33mpi for non-gauged defects. Although we claimed that the 
critical value of rj does not depend on initial conditions [13], it seems to lack of convincible 
explanations according to recent arguments [5-7]. We therefore return to the same model 
and clarify the effect of initial inhomogeneity both analytically and numerically. 

The plan of this paper is as follows. In Sec II, we discuss analytically the scalar held 
dynamics under the approximation that the metric is homogeneous. In Sec HI, we numeri- 
cally solve the Einstein equations for a global monopole with a large gradient. Section IV is 
devoted to a summary and discussion. We use the units of c = ^ = 1 throughout the paper. 

II. SCALAR FIELD DYNAMICS IN A HOMOGENEOUS SPACETIME 

Let us consider the scalar fleld dynamics under the approximation that the metric is 
homogeneous, following Berera and Gordon [7]. The equation of the scalar fleld ^(i,x) is 



where an overdot denotes time-derivatives, a{t) is the scale factor and H = a/a. Although 
a{t) is not a solution of the exact Einstein equations, it is instructive to write down its 
evolution equation under the present approximation. 
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This equation tells us that inflation does not take place (a < 0) if the kinetic term dominates 
over the potential. It should be noted that a gradient term does not appear in Eq.(2.2). 
Gradient energy can affect cosmic expansion only by exciting kinetic energy through the 
field equation (2.1). 

Because inflation can take place at 0, we approximate the potential (1.1) as V — 
\r]'^/A - \rf'^'^/2. By Fourier expansion ^'(t,x) = E *k(^)e*'', Eq.(2.1) becomes 

^k + Sii'^k+^^k-ATy^^k^O. (2.3) 

We assume the slow-roll condition for a homogeneous and isotropic spacetime, |^'| ^ 3i/|\l'| 
and ^^/2 ^ y at ^ 0, which is equivalent to 77 » mpi/Vdn. (Actually the numerical 
flnding r) > 0.33mpi [13] is more precise condition.) Under this condition what happens 
if the A:^-term of (2.3) dominates? Clearly ^k begins oscillation, which contributes kinetic 
energy. In the case of other nontopological inflationary models, inflation never happens if 
kinetic energy dominates until the mean field rolls down to the potential minimum and the 
slow-roll condition breaks down. On the other hand, in the case of topological infiation, the 
false vacuum ^ = is topologically stable and never decays. As long as the core region 
\& is not contracting at an initial time (d > 0), oscillations are eventually damped. 
Because a local region of ^ ^ always exists, inflation takes place eventually. 

The above arguments are rough because we have neglected the inhomogeneity of the 
spacetime. This is why we shall move on to numerical analysis in the next section. Never- 
theless, the above arguments will help us to understand our numerical results intuitively. 



III. NUMERICAL ANALYSIS OF A MONOPOLE WITH A LARGE GRADIENT 

Here we take inhomogeneity of spacetime into account by solving the Einstein equations 
coupled to the scalar held equation numerically. We restrict ourselves to the global monopole 
system. We believe that the main features of dynamics do not change for domain walls and 
global strings, according to the previous work [13]. Because we have already known that 
the inhomogeneity of gauge flelds affects the onset of inflation [14] and our present concern 
is pure effects of the scalar fleld inhomogeneity, we do not consider gauged defects here. 

We therefore consider the Einstein-Higgs system, which is described by the action. 



^ — J d^x^f^ 



^7e--(M»)2-i^(^) 



1,2,3), 



(3.1) 



with the potential (1.1). The action (3.1) yields the field equations. 



2 mpi 



(3.2) 
(3.3) 



with 



7;, = V^^^V.*'^ - g^. 



(3.4) 
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We assume a spherically symmetric spacetime and adopt the metric, 

ds" = -de + A\t, r)dr'' + rYide'' + sin" ^V)- (3-5) 
For the scalar field, we adopt the hedgehog ansatz: 

= ^(t, r)(sin 6 cos ip, sin ^ sin 99, cos 9) (3.6) 

In the following we shall write down the field equations (3.2) and (3.3) with the assump- 
tions (3.5) and (3.6). To begin with, we introduce the extrinsic curvature tensor Kij-. 



A 



K = Kl 



Next, following Nakamura et al. [15], we introduce auxiliary dynamical variables. 



A-B 



a = 



k = 



K'l - k: 



(3.7) 



(3.8) 



and use a new space variable, x = r^, instead of r. The advantages of those variables were 
discussed in Ref. [16]. 

We introduce further auxiliary variables, 



B' 



i) = —, w = ijj, ^ = ip', 



(3.9) 



where a prime denotes d/dx. A full set of dynamical variables is a, S, C, tp, ^, k and 
VD. The time-derivatives of the first five variables are given by the definitions above: 

d^ -aKl + Bk, B^-BK^, C ^ Kf , (3.10) 

iP = w, i^w'. (3.11) 

Note that the value of K^' in Eq.(3.10) is determined by the momentum constraint (3.13) 
below, which gives a more accurate value than a finite difference of Kq. 
Now we write down the field equations (3.2) and (3.3) as 
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k-kK + 
Stt 



nipi^A^ 



A^ 

A^{x^ + V') 



2AC' - 2A'C + aC + a' 
ilj^a{A + By 



(3.12) 
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(3.14) 
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(3.15) 



zu = Kw + 

2ij 
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Equations (3.12) and (3.13) are constraint equations, which fixes initial values of the metric. 
Equations (3.14)-(3.16) together with (3.10) and (3.11) provides the time evolution of the 
seven variables. 

Let us discuss how to discretize a space at each time into a mesh and approximate spatial 
derivatives. Because our field equations are written with x instead of r, the simplest way is 
to use a regular mesh with separation Ax in a;-space, as we did in Ref. [16]. However, this 
mesh gives lower resolution for a center than for a far region in terms of r-space, which is 
unsuitable for the present case, where the central region has a large gradient. Therefore, we 
adopt a regular mesh with separation Ar in r-space, 

n = 0, r2 = Ar, n^{i- l)Ar, (3.17) 

and approximate x-derivatives of any function F{x) by 

dF 1 dF 1 Fj+i - Fj+i - /o io\ 

—^{Xi) — l-{Xi) — ; = . (3.18 

dx ^ ' 2ri dr ^ ' 2ri 2Ar x^+i - Xi_i ^ ' 

We impose further regularity condition on all dynamical variables at r = 0, F(0), as 
follows. Although all evolution equations do not contain any diverging factor at r = like 
1/r, we do not use them to obtain Fi = F(0) at each time. Instead, following Hayley and 
Choptuik [17], we employ a "quadratic fit", 

4-F2 - Fg , , 

Fi = '3 (3.19) 

We find that this simple method is very effective to keep good accuracy in the center. 
As initial conditions of the scalar field, we suppose 

m(t = Q,r) = r] tanh , ^(t = 0, r) = 0, (3.20) 

where 5 is a free parameter, which controls the initial gradient at ~ 0. If the previous 
arguments applied to topological inflation, 5 would have to be larger than the horizon size 
Hq ^ = (87r\/(0)/3mpi2)-i/2 for the onset of inflation. Now our question is specifled: can 
inflation start even ii 6 <^ Hq^7 We therefore choose 6 = O.OIHq^ below as an example. 

For the metric, there are four unknown variables, B, a, K and A;, while there are two 
constraint equations (3.12) and (3.13). Therefore, two of the four variables are arbitrary 
chosen, and this choice is an important problem. If we choose K =constant, for example, 
large initial inhomogeneity of ^ causes a black hole or a singularity, as discussed by Chiba 
et al. [3] The appearance of black holes or singularities, however, does not mean failure 
of inflation, because inflation is a local phenomenon. To judge whether inflation occurs or 
not, we should solve the evolution equations; however, it is difficult to solve them for such 
spacetimes numerically. For such a technical reason, we set initial conditions without a 
black hole nor a singularity by adopting the following two ways, under the condition that 
the initial 3-space is not contracting, < 0. 

(A) Following Goldwirth and Piran [1], we set a homogeneous and isotropic 3-space by 
introducing another massless scalar fleld x to compensate inhomogeneous energy density of 
^. We assume x(^ = 0,r) = and x{t = 0,r) > 0. The initial metric functions become 
B — 1, a = 0, k — Q and K ^constant. 
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(B) WithoTit introducing a massless scalar field, we set B — 1 and a — and solve the 
constraint equations (3.12) and (3.13) to obtain K and k. 

We demonstrate solutions with A = 0.1, r] = mpi and 5 = 0.01i7(7^ in Figs. 1 and 2. 
In both cases we see in (a) and (b) the monopole core expands. As Eq. (3.14) suggests, 
whether the expansion rate is increasing or not is determined by the sign oi "if^ — V even 
in an inhomogcncous spacctime. Therefore, once kinetic energy becomes smaller than the 
potential, inflation starts even though gradient energy is still dominant. This behavior is 
actually seen in Figs. 1(c)- (e) and 2(c)-(d). 

The evolution time in Fig. 2 is relatively short because calculation breaks down soon 
after the black-hole horizon appear at t ^ 0.27Hq^. However, we beheve it demonstrates 
sufficiently that inflation actually begins. 

It should be emphasized that topological inflation can take place even if Pgrad ^ V, as 
shown in Figs. 1(c) and 2(c). This is in contrast with the previous results [1], which showed 
that new inflation requires a more stringent condition than Pgrad < V- 

IV. SUMMARY AND DISCUSSION 

We have investigated whether topological inflation requires initial homogeneity over the 
horizon size. Contrary to the results for nontopological inflationary models, which have 
been investigated by many people, topological inflation can take place even if the initial 
false vacuum has one percent of the horizon size and gradient energy of 10^V"(0). This is 
due to topological stability of false vacuum. In this sense topological inflation explains most 
naturally the onset of inflation in a really chaotic universe. 

Here we make a comment on the papers by Vachaspati and Trodden [5], who considered 
initial conditions for inflation based on the null Raychaudhuri equation, and concluded that 
homogeneity on super-horizon scales must be assumed as an initial condition. They claimed 
that the size of the initial inflationary patch, X, must be greater than the inflationary 
horizon, -f^j^j, which must be larger than the (pre-inflationary) background FRW inverse 
Hubble size at the time inflation starts, Hp^^^, that is, 

X > Hr) > i/-!^. (4.1) 

However, what they really proved is only the second inequality, H~^jr > Hp]^-^, to which we 
have no objection. On the other, they just assumes the first inequality, although its validity 
is more crucial. Actually, what we showed in the present paper is topological inflation can 
takes place even \i X < 

Another note is added. Stabihty of false vacuum is a key feature of topological inflation, 
which sometimes leads to the idea that inflation cannot end [12]. It is, however, not to worry. 
Once inflation begins, the evolution of a local region is described by Eq.(2.1) without the 
gradient term. Therefore, just like new inflation or chaotic inflation, ^' simply rolls down 
to the potential minimum except for = 0, that is, inflation ends in any region except the 
one point. More detailed arguments are given in Ref. [18]. 
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FIGURES 
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FIG. 1. A solution for the case (A). We show the evolution of the inflaton field, the volume element, 
and gradient energy Pgrad = (3\l'/5r)^/(2A^) in (a), (b) and (c), respectively. For kinetic energy 
Pkin = {d^ /dt)'^ /2, (d) and (e) report its early and late behaviors, respectively. The abscissa is 
the arial radius of the spacetime B{t,r)r, normalized by the horizon size Hq^. 
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FIG. 2. A solution for the case (B). We show the evolution of the inflaton field, the volume 
element, gradient energy and kinetic energy in (a), (b), (c) and (d), respectively. 
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